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Abstract 

The negatively-charged nitrogen-vacancy (NV ) center in diamond is at the frontier of 

quantum nano-metrology and bio-sensing. Recent attention has focused on the application 

of high-sensitivity thermometry using the spin resonances of NV - centers in nano-diamond to 

sub-cellular biological and biomedical research. Here, we report a comprehensive investigation 
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of the thermal properties of the center’s spin resonances and demonstrate an alternate all- 
optical NV thermometry technique that exploits the temperature dependence of the center’s 
optical Debye-Waller factor. 

The negatively-charged nitrogen-vacancy (NV - ) center in diamond- 1 is an important defect 
system for a range of quantum technologies. The center has been key to several advances in quan¬ 
tum information processing, including the realization of solid-state spin quantum registers and ele¬ 
ments of quantum networks The center has also attracted significant attention due its potential 
as a high-sensitivity nanoscale quantum sensor that can operate in ambient and extreme condi¬ 
tions. Indeed, the NV - center has been employed in many impressive demonstrations of nano- 
magnetometry, 7-10 electrometry,----’ 12 piezometry 13 and thermometry.— — 8 Furthermore, the cen¬ 
ter may be incorporated into nano-diamonds, whose chemical inertness and biocompatibility, allow 
sensing techniques to be performed within living cells. 19 Recently, the center’s bio-sensing appli¬ 
cations were expanded with the demonstration of sub-cellular temperature gradient mapping and 
control using NV - nano-diamond and gold nano-particles within a human embryonic fibroblast. 17 
This exciting demonstration can be potentially extended to in vivo thermometry and thermoablative 
therapy.— 

Each of the center’s metrology applications either directly employ its electron spin or a nuclear 
spin that is coupled to its electron spin as a quantum sensor. High sensitivity is principally achieved 
by the unique combination of the exceptional electromagnetic, mechanical and thermal properties 
of diamond and the long-lived coherence of the center’s electron spin, which persists in ambi¬ 
ent and extreme conditions.- 14 *^ 0,21 Nanoscale sensing is enabled by the atomic size of the NV - 
center, its bright fluorescence that allows single centers to be located with nano-resolution and its 
mechanism of optical spin-polarization/ readout (OSPR) that allows the magnetic resonances of its 
electron spin to be optically detected (ODMR).- 

It is evident that a thorough understanding of the effects of temperature on the NV - center’s 
electron spin resonances is required to optimize its implementation in quantum technologies for 
ambient/ fluctuating environments.— It is also clear that specific understanding of these effects in 
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nano-diamond NV centers (as opposed to bulk diamond) is required for the refined realization of 


the center’s bio-sensing applications. Such specific differences of nano-diamond include variations 
in thermal properties and generally higher intrinsic strain and abundance of local charge traps. 
Despite these requirements, only the temperature shift of the zero-field (no strain) ground state 


14.23-26 


electron spin resonance has been examined in detail. 


Consequently, the center’s excited 


state spin resonances represent additional quantum resources for the center’s exciting applications 
that are yet to be precisely examined. 27 

To rectify this situation, we report comprehensive observations of the effects of temperature 
on all of the optically detected spin resonances of NV “ centers in nano-diamond. Using models 
derived from first principles, we describe the temperature variations of each of the zero-field and 
strain interactions that define the electron spin resonances. In doing so, we provide significant 
insight into the thermal properties of the resonances and the influence of the nano-diamond envi¬ 
ronment. Additionally, we demonstrate an alternate all-optical NV~ thermometry technique that 
exploits the temperature dependence of the center’s optical Debye-Waller factor. Our all-optical 
technique has a projected temperature noise floor of 0.1 K Hz -1 / 2 for a nano-diamond containing 
500 NV~ centers. This noise floor potentially surpasses that of other existing biocompatible nano¬ 
thermometry techniques, except for the NV~ ODMR technique. I7 ' 28 70 Indeed, the simplicity and 
robustness of our all-optical technique may prove advantageous compared to the ODMR technique 
for particular biological applications where microwave excitation is prohibitive. 

The NV~ center is a Cj v point defect in diamond consisting of a substitutional nitrogen atom 
adjacent to a carbon vacancy that has trapped an additional electron (refer to Figure |T^). Figure 
[lb depicts the center’s electronic structure, including the low-temperature zero phonon line (ZPL) 
energies of the visible (E\/ rsj 1.946 cV)~ and infrared (T)r ~1.19 eV}22 24 transitions. The spin 
resonances of the 3 A 2 and 3 E levels can be each optically detected at room-temperature as a change 
in the visible fluorescence intensity and described by a spin-Hamiltonian of the form 


H 



( 1 ) 
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where S and I are the dimensionless spin-1 operators of the electron and 14 N nuclear spins, re¬ 
spectively, that are defined with respect to the coordinate system depicted in Figure [lji, and the 
parameters {Q), A, srf 3 -) differ between levels and are identified in the following text. 

As depicted in the inset of Figure [TJd, the ground 3 A 2 level exhibits a zero-field fine structure 
splitting between the m s — 0 and ±1 spin sub-levels of ^ g . s . ~ 2.87 GHz (room-temperature) due 
mainly to electron spin-spin interaction.— Under crystal strain that distorts the C 3 V symmetry of 
the center, the m s = ±1 sub-levels are mixed and their degeneracy is lifted . 16 This strain dependent 
splitting of the m s — ± 1 sub-levels is 2Ag. s .. The 3 A 2 level also exhibits 14 N magnetic hyperfine 
structure described by the parameters ' e ^g!s. —2.14 MHz and ~ —2.70 MHz, 37 but this 

minor additional structure does not feature in this work. Nor do the 14 N electric hyperfine inter¬ 
actions. At low temperatures (< 10 K), the excited 3 E level exhibits a non-trivial fine structure 
arising from a combination of electron spin-orbit, spin-spin and strain interactions. 38 Above 150 
K, 38 the observed fine structure of the 3 E level becomes analogous to that of the ground 3 A 2 level, 
with a zero-field fine structure splitting between the m s — 0 and ±1 spin sub-levels of f^ e .s. ~ 1-42 
GHz (room-temperature) and a strain dependent splitting 2<f es . of the m s — ± 1 sub-levels (see Fig¬ 
ure [B- 1940 Above 150 K, the 3 E level also exhibits 14 N hyperfine structure analogous to that of 
the 3 A 2 level, but with much larger parameters rs./ ^g.s. rs./ 40 MHz (assumed isotropic)— and 
is consequently important to this work. The collapse of the observed 3 E fine structure above 150 K 
is thought to be due to phonon mediated orbital averaging, however a precise model of this process 
is yet to be presented.— Consequently, the physical origins of the strain dependent <%. s . parameter 
have not been previously established. 1 The 3 E hyperfine structure has also not been previously 
examined theoretically and it has not been observed at low temperatures. - 

In our experiments, we measured the temperature dependence of the ODMR of NV centers 
in type lb HPHT nano-diamond crystals (see Figure [I])- On average, the nano-diamonds had a 
diameter of ~ 30 nm and contained ~ 15 NV~ centers. We performed ODMR measurements on 
a total of 10 nano-diamonds. Photoluminescence (PL) spectroscopy was also performed on five of 
these nano-diamonds and the results of these five nano-diamonds were analysed and found to be 
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Figure 1: (a) Schematic of the NV center depicting the vacancy, the substitutional nitrogen atom, 
the surrounding carbon atoms and the adopted coordinate system (z axis aligned with the C^ v axis 
of the center and the x axis is contained in one of the center’s mirror planes), (b) Schematic of 
the center’s electronic structure, including the low-temperature visible Ey ~ 1.946 eV and infrared 
£ir ~1.19 eV ZPL energies. Lower box - The fine structure of the 3 A 2 level: at zero field with a 
single splitting @> g s . ~ 2.87 GHz (room-temperature); and under symmetry lowering strain, with 
an additional strain dependent splitting 2<f gs .. Upper box - The fine and hyperfine structures of 
the 3 E level observed at room-temperature: at zero field with fine ^ e .s. ~ 1-42 GHz and hyperfine 
~ 40 MHz splittings; and under strain, with the additional splitting parameter Ae. s .. (c) the 
experimental setup depicted in two sections: (top) the front projection depicting the laser spot 
partially covering the U-shaped gold wire connected to a RF generator; (bottom) optical diagram. 
Labels are defined as: O - microscope objective, D - photodetector, L - laser, F - optical filter and 
S - sample. 
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consistent. The ODMR results for one representative nano-diamond are presented in the following, 
whereas PL results for several nano-diamonds are presented as specified. The nano-diamonds were 
spin coated on a fused quartz substrate. The spin resonances were driven by microwaves carried by 
a U-shaped, 4 cm long and 200 x 0.5 /dm in cross-section band of deposited gold on the substrate. 
The investigated crystals were located at distances between 3 and 30 jUm from the edge of the gold 
band. The spin resonances were optically detected using 532 nm laser excitation and fluorescence 
collection via an epifluorescence design. ODMR spectra were measured using continuous wave 
(c.w.) optical and microwave excitation. As the microwave field was weak and the crystals had 
good thermal contact with the substrate, the crystals were primarily heated (up to 600 K) by the 
laser light illumination of the gold band. See supplementary information for further experimental 
details. 

Our all-optical nano-thermometry technique uses the Debye-Waller factor (DWF) of the NV 
visible transition. The DWF is the ratio of the area under the ZPL and the total emission band (see 
Figure 12). The technique was calibrated by measuring the DWF of nano-diamond NV~ centers 
inside an oven operating under stabilized temperature conditions (see Figure [ 2 - For T <C 7b, the 
temperature dependence of the DWF is well described by the simple Debye model of the phonon 

continuum— 


DWF = (2) 

where 7b is the Debye temperature of the host crystal (7b ~ 2220 K for bulk diamond) and S is 
a parameter defining the electron-phonon coupling strength. For the nano-diamonds in the oven, 
the observed value of 7b was ~ 1614(23) K, which is 1.37 times smaller than the value for bulk 
diamond derived from its thermodynamic properties at low temperatures. Such a mismatch is typ¬ 
ical and is well established in other materials.— Given the observed value of 7b, the temperature 
of the nano-diamonds subjected to the laser-assisted heating was found to follow a linear function 
T = To + bP\ as of laser power P las , where Tq ~ 294 K is the room-temperature and the proportion- 
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ality constant b was found to be 0.51(3) K/ mW for the crystal whose data is shown in Figure [2]:. 
The electron-phonon coupling parameter 5 was determined separately for the nano-diamonds in 
the oven and those used for the ODMR experiments. 



Wavelength (nm) 

(a) 



Temperature squared f f /K ? ) x i q 6 

(b) 



(c) 


Figure 2: Example NV visible ZPL spectrum (a) at room-temperature and the DWF as a function 
of (b) oven temperature T and (c) laser power F\. dS . (a) demonstrates the fit (red solid) of the ZPL 
(blue points) using two strain-split lorentzian lineshapes (black solid) and a linear function (red 
dashed) to account for the contribution of the phonon sideband underneath the ZPL. The equation 
of the straight line in (b) is logDWL = —5(1 + |t 2 T 2 /Tfi), where 5 = 4.57(7) and To = 1614(23) 
K are derived from a least-squares fit to the data points (circles). The dependence of the DWL 
on the laser power in (c) is described by the substitution T = Tq + bP\ dS ,, where Tq — 294 K is 
the room-temperature and a least-squares fit of the data points (circles) yields 5 = 4.79(6) and 
b = 0.51(3) K/ mW. In (b) and (c), error bars are smaller than the point size (see supplementary 
information). 


Given the collection of N photons emitted by NV~ centers in a nano-diamond, the smallest 
temperature change that can be detected by the optical DWL thermometry technique is defined by 
the shot-noise uncertainty in the measurement and the temperature dependence of the DWL 


STrr 


DWL -|- v'DWF /c/DWF \ 

V dT ) 


v'DWF /JDWL\ 1 
VN \ dT ) 


(3) 


The presence of background fluorescence raises the uncertainty by a factor of y/\ +3r, where 
r is the ratio of the approximately uniform background intensity under the NV~ ZPL and the 
ZPL peak intensity. 46 The value of N is related to the measurement time T, the optical collection 
efficiency /i, the single NV - photon emission rate y and the number of NV centers n via N — 
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nfj,yx. Consequently, the temperature noise floor r\j of the DWF thermometry technique is 


Vr = 8T min Vr = \/(l +3r)C Z p/ 2 fl>, (4) 

where Czpl = n/i /DWF is theZPL emission rate and <t> = DWF(r/DWF/c/7) At room-temperature, 

exhibited little variation over our nano-diamond sample and had an average value of <t> = 
154(10) K. Using the typical optical lifetime of NV centers in nano-diamond— and the details 
of the optical setup,— we can estimate y ~ 40 MHz, /i ~ 0.021 and DWF~ 0.005. The estimated 
noise floor of a single NV “ center at room-temperature and in the absence of background fluores¬ 
cence is thus T]t ~ 2. 3 K Hz” 1 / 2 . If n 500 as in the recent bio-thermometry demonstration,— 
the the noise floor becomes r/r ~ 0.1 K Hz ! / 2 . 

To test our estimate of the noise floor of the DWF thermometry technique, we characterized 
the noise in our DWF measurements and performed a series of temperature measurements. Figure 
[3jr demonstrates the Poissonian fluctuations of a fluorescence spectrum of a nano-diamond. This 
spectrum was collected over one second and contains fluorescence from NV - centers as well some 
neutral NV (NV°) centers present in this particular nano-diamond. The upper panel of figure [3b 
depicts a sequence of DWF measurements, including a step caused by a 17 K temperature increase 
via laser heating. Each measurement was averaged over one second. The observed standard de¬ 
viation of the DWF measurements corresponds to a temperature uncertainty of 4.0(4) K, which 
directly implies a temperature noise floor of 4.0(4) K Hz -1 / 2 . The lower panel of figure [3b depicts 
a series of temperature measurements using a second nano-diamond at room-temperature over a 
period of 23 minutes. This nano-diamond contained more NV - and NV° centers and each mea¬ 
surement was averaged over 50 s. The standard deviation of the measurements from the mean value 
was 0.50(5) K, which corresponds to a noise floor of 3.5(3) K Hz '/ 2 . The deviation from a cubic 
polynomial trend of the room temperature is even smaller, 0.3 K. The noise floors achieved using 
these two nano-diamonds are in agreement with (4). This can be confirmed using the experimen¬ 
tally observed values of Czpl ~ 14 kHz and 42 kHz and r ~ 3.7 and 10.5 for these nano-diamonds, 
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respectively. As seen in figure [3ji, the non-zero NV° fluorescence makes the main contribution to 
the background. The increased noise floor arising from the NV° background can be reduced in 
several ways. Since the NV° center ZPL occurs at 575 nm, NV° fluorescence can be avoided if 
optical excitation is instead performed at wavelengths in the range 575 — 637 nm. Alternatively, 
nano-diamonds with relatively few NV° centers could be purposefully fabricated or post-selected. 
Hence, NV° background fluorescence does not present a fundamental barrier to achieving the pro¬ 
jected noise floor r\j ~ 0.1 K Hz -1 / 2 for a nano-diamond containing ~500 NV~ centers. 

Wavelength (nm) x 10 ' 3 





(a) (b) 

Figure 3: (a) Upper panel: PL spectrum collected over one second at room-temperature (up¬ 
per panel) of a nano-diamond containing NV and NV° centers. The ZPLs of the NV charge 
states are as denoted and the NV~ ZPL area is outlined in red. Lower panel: the standard 
deviations of the number of counts per spectrum pixel Nj from its mean value M, determined 
by a spectrum collected over a much longer time (50 s). The equality of the standard devi¬ 
ation of (Nj — Mi) / y/M~i to 1 demonstrates the Poissonian noise of the spectrum, (b) Upper 
panel: Time series of DWF* measurements that were each obtained after one second of aver¬ 
aging. DWF* is distinguished from DWF and is defined as the NV~ ZPL area divided by the 
area of the total fluorescence spectrum (including the background NV° fluorescence). Note that 
DWF*(r/DWF* /dT)~ l = DWF(dDWF/dT)~ l = <£>. A temperature step of 17(2) K above room- 
temperature (as measured via DWF technique), achieved via laser heating, is evident at t — 23 
s. Lower panel: Time series of DWF temperature measurements (points) over 23 minutes under 
ambient conditions with a time-averaged temperature (dashed line) of 294 K and 0.5 K fluctua¬ 
tions. The solid cubic-polynomial line indicates possible long-time trend. Each measurement was 
averaged over 50 s. 


This projected sub-Kelvin noise floor of the DWF thermometry technique is comparable to 
the demonstrated 130 mK Hz 1//2 noise floor of the ODMR technique in nano-diamond— and it 
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potentially surpasses the noise floors of other biocompatible nano-thermometry techniques with 
sensor sizes less than 100 nm. ——The notable competitor being CdX (X= S, Se or Te) quan¬ 
tum dots, whose noise floors have not been clearly determined, but based upon some demon¬ 
strations, could rival both NV nano-thermometry techniques.■ 29,3 ° The projected noise floor of 
the ODMR technique for ideal nano-diamonds containing 500 NV~ centers is ~1 mK Hz -1 / 2 .— 
Whilst this is much lower than the noise floor of the DWF technique, the fabrication of such ideal 
nano-diamonds, with negligible strain inhomogeneity and sufficient purity to support NV - spin 
coherence times of ~ 1 ms, is yet to be achieved. Comparing further with the ODMR technique, 
the DWF technique is a much simpler and more flexible all-optical technique that does not require 
microwave control, and thus may be better suited to particular biological applications. Indeed, the 
DWF technique can be immediately implemented using a commercial spectral-resolving imaging 
system, such as a Raman microscope. Like the ODMR technique, the DWF technique has the po¬ 
tential to be extended to in vivo applications if a different method of fluorescence spectroscopy was 
employed, such as two-photon excitation spectroscopy. 49,50 Spatial resolution may also be further 
improved by all-optical far-field sub-diffraction imaging techniques. 51,52 Thus, there is clear moti¬ 
vation for further investigation of the DWF thermometry technique as an altemate/complementary 
technique to the ODMR technique and other current nano-thermometry techniques. 

There has been several studies of the effects of pressure and temperature on the visible ZPL, 
the infrared ZPL and the zero-field 3 A 2 spin resonance S> g , s , of centers in bulk diamond. 14 - 2 - 1 - 26 - 53 
Acosta et al have also studied the temperature dependence of the 3 A 2 strain splitting S{, s in bulk di¬ 
amond NV - ensembles, but since it was much smaller than did not observe a variation above 
experimental uncertainty over the range 5-400 K. 2 -' 24 - The pressure and temperature dependence of 
the 3 E splittings D s s and <f e . s have not been previously studied. Recently, Doherty et al 26 demon¬ 
strated that the model original developed by Davies 53 to explain the temperature shift of the visible 
ZPL energy, successfully describes the temperature shifts of the infrared ZPL energy and In 
Davies’ model, there are two origins of the temperature shifts of the spin resonances: (1) the spin 
energies are perturbed by the strain of thermal expansion, and (2) the vibrational frequencies as- 
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sociated with different spin states differ. The later is a consequence of quadratic electron-phonon 
interactions.- 26 Since thermal expansion and hydrostatic pressure are intimately related, the con¬ 
tribution of thermal expansion to the temperature shifts may be alternatively expressed in terms 
of the hydrostatic pressure shifts and the pressure of thermal expansion in diamond. 53 Thus, the 
temperature shift of the 3 A 2 zero-field resonance Af^ g . s .(7) = A^| x s (T) + AfZ|. s p (T) is the sum of 
thermal expansion ( T) and quadratic electron-phonon A^|' s p (7) contributions given by 

a = r g ,.P(T) 

A^|“ p (r) = jf n(co,T)8 g . s X(o)p{cj)d(0, (5) 

where r g . s . = 14.58(6) MHz/GPa— is the hydrostatic pressure shift of S> g . s ., P(T ) = B Jq e(t)dt 
is the pressure of thermal expansion, B = 442 GPa is the bulk modulus of diamond, e(T) is the 
diamond volume expansion coefficient, n ((0 , T) — (e l,0) / k,il — 1)~ 1 is the Bose-Einstein distribution 
of vibrational occupations, h is Planck’s constant, kg is Boltzmann’s constant, Cl is the highest 
vibrational frequency of diamond, p(co) is the vibrational density of states, and 5 giS . (to) is related 
to the differences in the vibrational frequencies of the 3 A 2 electron spin states. 26 Given power 
series expansions of e(T) and <5 g . s .(ft)jp(ft>), the temperature shift A 3> g . s .(T) may be approximated 
by a polynomial in temperature 7. 26 Doherty et al’s application of Davies’ model is extended 
here to ff g s . and ^ e s . by deriving analogous expressions for A S> e .s.(T) an d AS g s (T), where the 

corresponding T and S(co) parameters differ from A S> g & (T). Explicit expressions for each of the 

r and 8 (to) parameters can be obtained from first-principles via the application of the molecular 
orbital model of the center (see supplementary information). 54i5S 

A model of the collapse of the 3 E fine structure above 150 K and the temperature dependence 
of S’e.sX'P) can be derived from the low temperature 3 E fine structure in the condition where spin- 
conserving electron-phonon transitions within the 3 E level occur much faster than the frequencies 
of the 3 E fine structure interactions and the optical decay rate (see supplementary information). 
Observations of the temperature dependence of the 3 E ODMR intensity— and the visible ZPL 
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width— imply that this condition is satisfied at temperatures > 150 K. Recognising that the 3 E fine 
structure can be described by the interaction of orbital and spin sub-systems, this condition implies 
that thermal equilibrium of the ?1 E orbital states is achieved within the lifetime of the 3 E level and 
that no coherences exist between the orbital and spin sub-systems. 57 Hence, the orbital and spin 
sub-systems may be decoupled, which leads to the expression 


4.S XT) = Di,M(T) 


( 6 ) 


where D e ; s ~ 0.775 GHz is the low temperature transverse electron spin-spin interaction of the 3 E 
level— and M(T) = (e h ^' i ktiT — 1)/(V'? A/A + 1) is a strain-temperature reduction factor derived 
from the Boltzmann distribution of the orbital sub-system and is the strain splitting of the 3 E 
fine structure. The above expression clearly demonstrates that the physical origins of <f e . s . is a 
product of transverse spin-spin interactions D~ s and a factor 3% that explicitly depends on the ratio 
of the 3 E strain interaction and temperature. In addition to the explicit factor 3%, <f es . is intrinsically 
dependent on strain and temperature via their effect on the spin-spin interaction D^ s . This intrinsic 
dependence will be analogous to the previously described strain and temperature variations of the 
other parameters ^ e .s. and<f gs .. 

The observed temperature variation of the ground state 3 A 2 spin parameter , s . is depicted in 
Figure @Ji. Referring to Figure 0^, it can be seen that over the majority of the temperature range, 
the temperature shift A^ g s .(r) = £^ giS .(r) — ^ g . s .(294K) observed here in nano-diamond is very 
similar to the previous observation in bulk diamond, — thus implying that there is little variance 
of A^ g s (T) in nano-diamond. This conclusion is important to the large-scale implementation 
of NV~ nano-thermometers based upon A^ g s .(r) since it suggests that there will be little varia¬ 
tion between nano-thermometers. 2 - Furthermore, as temperature in our nano-diamond A S> g s (T) 
measurements was calibrated using the optical DWF thermometry technique, the similarity of our 
observations to those in bulk diamond further validate the optical DWF technique. 

As experienced by Acosta et al in the bulk diamond measurements over temperature 5-400 
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(a) (b) 

Figure 4: (a) The observed temperature variation of ^ g . s . and the fit of a quadratic polynomial 
^g.s. = dg.s.+bg.s.T + c g .s T 2 , where a g . s . = 2870(3) MHz, b g . s . = 6(1) x 10~ 2 MHz/K and c g . s . = 
—2.3(2) x 10 -4 MHz/K 2 . (b) Comparison of the temperature shift A 3> g . & (T) observed here in two 
different nano-diamonds (blue and red solid) and the previous observation in bulk diamond (black 
dashed).— Error bars in (a) are smaller than the point size (see supplementary information). 


K, 23-24 we did not conclusively observe a temperature variation of <f g . s . up to 600 K. The difficulty 
in observing the temperature shift of ^ g . s . arises from its much smaller magnitude (~ 10 MHz here), 
which based upon the 0.01 fractional change of ^ g s . over the range 300-600 K, implies a variation 
of <4g. s . at the limit of detection (~ 0.1 MHz here). The first-principles model of <C s . identifies 
that the temperature variation of 4 g . s . is potentially much more complicated than ^ g s . due to a 
mixture of strain and transverse spin-spin factors (see supplementary information for discussion). 
The approximate temperature independence of 4 gs in nano-diamond is, however, an important 
observation for ODMR thermometry techniques that utilize a single ground state spin resonance 
(with frequency Si ± A gs . in the absence of a magnetic field), since the temperature variation of 
the spin resonance is restricted to just that of f§? g . Si . 

At each temperature, the observed 3 E ODMR spectra exhibited just two lines with no hyperfine 
structure (see Figured). This is due to the strain at this particular center being sufficiently large 
that the hyperfine resonances overlap (as depicted in Figure [I])- In this regime, the minor splitting 
of the m s — ±1 sub-levels is the weighted average of the hyperfine resonances 2 e e .s.(T), where 


£e,.(n 


— S’g s 

3 


•< r > + 3 


2 r 


+ ( ?e.s.(T) 


(7) 


The observed temperature variation of f^ e .s. and £ e . s . are depicted in Figures [5]} and[5]:. As demon- 
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strated in Figure [5}i, the relative shift of is much less than S'g.s.- Furthermore, as seen in 
Figure Sx the temperature shift of ^ e .s. is well described by just the contribution A^ X (T) of 
thermal expansion with a hydrostatic pressure shift of r e . s . = 11(1) MHz/GPa. Given the observed 
3 E ODMR linewidths, this shift indicates a projected pressure sensitivity of ~ 8(1) MPa Hz -1 / 2 , 
which is comparable to the ~ 0.6 MPa Hz -1 / 2 sensitivity of the ground state ODMR — The ab¬ 
sence of contributions of electron-phonon interactions is theoretically interesting and also explains 
why the relative shift of 3> e is much smaller than 3>g. s .. The implications of these conclusions for 
the first-principles theory of the temperature shifts is explored in the supplementary information. 



(a) (b) 



(c) (d) 


Figure 5: (a) Example 3 E ODMR spectra at different temperatures, (b) The observed temperature 
variation of £^ e . s . and the fit of the thermal expansion contribution (T) = T e s P(T). (c) The 
observed temperature dependence of <f e . s . and the fit given by the expression (??) with the single 
fit parameter h£±. (d) Contrast of the relative temperature shifts of (solid) and fF e .s. (dashed), 
where relative shift is defined as [S> a {T) — S> a {294K)\ / S> a (294K) for a = g.s. or e.s.. In (b) 
and (c), the temperature error bars are smaller than the point size and the frequency errors bars 
are derived from the uncertainties of lorentzian lineshape fits of the ODMRs (see supplementary 
information). 


As seen in Figure [5fc, a good fit of the observed temperature variation of e e . s .(T) is obtained 
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using the expression (??). The only free parameter in the fit is the strain energy h^± — 4.7(3) meV. 
This parameter is also the splitting of the visible ZPL, which was measured independently via 
room-temperature optical spectroscopy to be ~ 4.3(2) meV. This small, one standard of deviation 
discrepancy between the independent measures, supports our model of the temperature dependence 
of the 3 E fine structure. The implication of the fit is that the temperature variation of <f e s . is primar¬ 
ily described by the strain-temperature reduction factor &(T), rather than the intrinsic temperature 
variation of the transverse spin-spin interaction . Application of first-principles theory derives 
many similarities in the expected temperature variations of Z)j] s and (see supplementary in¬ 
formation). Consequently, the insignificant variation of Dy s is consistent with the observation that 
the variation of Z)jj. s is much smaller than Z)| s and is predominately due to thermal expansion. 

This new knowledge of the fine, hyperfine and thermal properties of the excited state spin 
resonances provides the basis for a more comprehensive assessment of the resonances as addi¬ 
tional quantum resources for the center’s applications. For example, the excited state spin may be 
employed to perform fast, optically gated swap operations on a nuclear spin memory in ambient 
conditions. 27 Alternatively, the differences in the properties of the ground and excited state spin 
resonances may be exploited to achieve bi-modal sensing, where the spins are conditioned to sense 
different effects concurrently. For example, the pressure sensitivities of the spin resonances of both 
states are comparable, but the excited state spin resonances in the absence of stress are insensitive 
to temperature compared to the ground state spin resonances. Consequently, the superior tempera¬ 
ture sensitivity of the ground state spin could be combined with the comparable pressure sensitivity 
of the excited state spin to sense pressure and temperature by iterating pulsed-ODMR sequences 
on the spins in rapid succession. 
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